Weak localization and magnetoresistance in a two-leg ladder model 
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We analyze the weak localization correction to the conductivity of a spinless two-leg ladder model 
in the limit of strong dephasing <C r tr , paying particular attention to the presence of a magnetic 
field, which leads to an unconventional magnetoresistance behavior. We find that the magnetic field 
leads to three different effects: (i) negative magnetoresistance due to the regular weak localization 
correction (ii) effective decoupling of the two chains, leading to positive magnetoresistance and 
(iii) oscillations in the magnetoresistance originating from the nature of the low-energy collective 
excitations. All three effects can be observed depending on the parameter range, but it turns out 
that large magnetic fields always decouple the chains and thus lead to the curious effect of magnetic 
field enhanced localization. 

PACS numbers: 73.63.-b, 71.10.Pm, 73.20.Fz, 72.15.Gd 



I. INTRODUCTION 

The interplay of interactions and disorder in elec- 
tronic systems confined to move in one spatial dimen- 
sion has remained a fascinating theoretical problem for 
many decades. It has long been known that any electron- 
electron interactions in a clean one-dimensional wire de- 
stroy the Fermi-liquid which is paradigmatic of higher 
dimensions and transform the system to a Luttinger liq- 
uid which has no coherent fermionic modes .^1 On the 
other hand, an arbitrarily small disorder potential in the 
non-interacting system leads to the localization of all 
states^K^ leading to strictly zero conductivity at any tem- 
perature. Since these pioneering works, a lot of progress 
has been made in understanding the effects of both dis- 
order and interactions in one-dimensional systems (in 
the following, we shall restrict the discussion to repul- 
sive interactions). Early worksP^ focused on the en- 
hancement of the backscattering off impurities induced 
by the Luttinger liquid, an effect seen already with a 
single impurity]^ while recently technical advances have 
allowed the study of strong localization in the presence of 
interaction ;] 10 ! 11 ! - these works concentrate on quasi-one- 
dimensional structures, but the same ideas should apply 
to true one-dimensional systems. It has even been possi- 
ble to take the more traditional ideas of weak localization 
and electron dephasing^ and apply them success fully to 
the peculiarities of the one-dimensional geometry! 13 ! 14 ! 

One of the nicest experimental manifestations of weak 
localization is in the magnetoresistance; in the traditional 
setting the applied magnetic flux breaks the time-reversal 
symmetry and thus destroys the weak localization correc- 
tion to the resistance. A weak magnetic field meanwhile 
has relatively little effect on the Drude (diffusive) compo- 
nent on the resistance, and thus allows the quantum weak 
localization correction to be isolated and studied experi- 



mentally. This unfortunately can not be seen in strictly 
one-dimensional nanowires however, as in this case the 
magnetic-field may be gauged out completely and will 
have no effect (the Zeeman coupling between the mag- 
netic field and the spin of the electrons produces a rather 
different effeclf^J). However, if one turns from single- 
chain nanowires to double-chain nanostructures, the so- 
called ladder models, then interesting orbital effects of 
the magnetic field may be restored. 

While the two-leg ladder is often introduced theoreti- 
cally as a beautiful model intermediate between one- and 
two-dimensional systems, it is also rather ubiquitous in 
nature. For example, carbon nanotubes which experi- 
mentally show power-law scaling behaviour in conduc- 
tance typical of the Luttinger liquid 16 are theoretically 
described in the low-energy regime by a model equiva- 
lent to that of a two-leg ladder.^ Curious magnetoresis- 
tance effects have been seen in nanotubes^ - however 
here one should be careful, as the not-completely-trivial 
mapping from the original nanotube to the two-leg lad- 
der means that the external magnetic flux will couple to 
the low-energy theory in a more complicated way than 
for the simple planar two-leg ladder. There are more 
direct experimental manifestations of the two-leg ladder 
however - including semiconducting double nanowires,^ 
as well as materials which are structurally made up of 
such ladders, such as PrBa 2 Cu 4 O g where control of dis- 
order and measurement of magnetoresistance are made 
in experiments.^ One can even artificially create such 
structures using cold atoms in an optical lattice,^ the 
disorder being added as a speckle pattern,^! and the ef- 
fect of a magnetic field being seen on the neutral atoms 
via the effect of an artificial gauge field P3 

Quite generically, ladder models have relevant 
backscattering terms (<?i terms in the parlance of Lut- 
tinger liquid physics), meaning that even the clean sys- 
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FIG. 1. Two-chain ladder model, where dots denote atom 
sites and solid lines indicate possible hopping directions. 




tem has non-trivial correlations ! 24 l 25 l At low tempera- 
tures, this leads to a rather surprising response to even a 
single impuritypSl and many different possibilities in the 
presence of a disorder potential.^^U There has also been 
a lot of interest in magnetic-field induced phase transi- 
tions in the clean ladder systems at low temperatureP^ES] 
While these works set the scene for the present study, 
we will be interested in quite a different regime where 
we have temperature sufficiently high that we can a) ne- 
glect the gaps opened by the interaction backscattering 
<7i terms; and b) consider sufficiently strong dephasing 
that the disorder potential may be treated peturbatively. 
We will show that even under these conditions, the inter- 
play between the interactions, disorder, and low dimen- 
sions gives us curious and counter-intuitive results, such 
as non-monotonic magnetoresistance, and magnetic-field 
enhanced localization. 

The paper is organized as follows: in Section II, we 
will introduce the ladder model that we will analyze. In 
Section III we will then study the effect of weak local- 
ization and magnetoresistance in two ways: one via a 
phenomenological introduction of a dephasing time, and 
the second from a true interacting microscopic model. 
We compare the two methods, and show that the phe- 
nomenological approach captures most of the essential 
physics. The exception is oscillations in the weak local- 
ization correction, which have their origin in correlation 
effects in the ladder model that are equivalent to spin- 
charge separation. In Section IV, having calculated the 
magnetoresistance, we analyze and explain its properties. 
Some technical details are relegated to the appendices. 
Throughout the paper, we use units where h = 1. 



II. THE LADDER MODEL 

We consider a two-chain ladder of spinless fermions as 
depicted in Fig. [II with lattice spacing a x and a y . The 
ladder is pierced Dy a magnetic field in the z-direction, 
and the electrons feel an interaction and a disorder poten- 
tial. The total Hamiltonian we will consider is therefore 
a sum of three parts: 

H = Hkin + Hi n t j s + #dis- (2.1) 

We now consider the three terms separately. 





FIG. 2. Band spectrum of the linearized model. The + and 
— bands exhibit the same Fermi velocity, but have different 
Fermi momenta, split by Afc = fc^,- — - 



A. Kinetic term and magnetic field 

We take the continuum limit in the x direction, which 
amounts to linearizing the spectrum; the inter-chain hop- 
ping term t± then serves to split the bands (assumed both 
to have the same Fermi- velocity vf)- 



H kiD = I dx 



ft 1 1 



(2.2) 



Here, the ^ u n ( x ) are the electron operators on the two 
chains, n = 1, 2 denotes the chain index and /i = R, L = 
+, — indicates the chirality (right and left movers). 
The Hamiltonian is diagonalized by the transformation 
= + (7 * i u,2)/'V / 2 with a = ±, which yields the 

dispersion relations 



s*(fc) = v F k + at x , 



;(k) = -(v F k + at x ). (2.3) 



The band structure is plotted in Fig. [2j note that the 
only quantitiy which discriminates the a = ±-bands is 
their respective Fermi momentum k F ,±- 

The magnetic field is included in Eq. (2.2 1 via the 
vector potential by means of the minimal substitution 
dx d x + ieA x in the longitudinal direction; and the 
Peierls substitution tx —> tj_e~ le J Adl in the transverse 
direction, where the integral is performed over a rung. 
We apply the magnetic field perpendicular to the ladder, 
B = Be z and choose the gauge Ai^ = T(Ba y /2)e x for 
the two chains 1 and 2. After Fourier transforming to 
momentum space, the Hamiltonian is given by 



H 



kin 



dk 

2^ 



+ (fc-*/)*U*)*/*,2(aO) 

+ t x (<i(fe)*M,2W + /i.c.)]. (2.4) 
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We have introduced the magnetic flux density 
/ = Ba y /2Q , 



(2.5) 



with the magnetic flux quantum $o — h/2e. Note that / 
has dimension [length" 1 ], whereas flux is conventionally 
dimensionless. The reason for this unusual dimension is 
the linearization of the spectrum, which is equivalent to 
omitting the lattice spacing in the x-direction. Neverthe- 
less, we take / to be a convenient parameterization of the 
magnetic field and for convenience will call it magnetic 
flux. 



The Hamiltonian (2.4 1 is diagonalized by the transfor- 
mation 



where 



(2.6) 



(2.7) 



and b € [0, 1] is a dimensionless measure of the strength 
of the external magnetic field 



v F TTf 



v/(W) 2 + *i' 



(2.8) 



This transformation leads to similar dispersion relations 
as in Eq. JO): 



4(k) = v F k ± y/(v F nf) 2 +tl, 



(2.9) 



The magnetic field affects only the Fermi momenta but 
leaves the Fermi velocity unchanged. The Fermi mo- 
menta splitting changes as 



k, 



kp + = Ak = 



i'F 



yJ(v F nf)*+tl. (2.10) 



Using Eq. (2.4 1, we also find a mapping between the 



Green functions in the chain basis and those in the band 
basis: 



G (i,i) 

(7(2,2) 

fj, 

G (l,2) 

(J, 



= ^[(l-/x&)G„, + + (l + M 6)G M ,_]. 
= i[(l + //&)G„, + + (l- M &)G M ,_]. 



G,, 



(2.11) 



Here, 



GM\tf;xa/) = -i(T t ^^{t,x)^{t',x')) (2.12) 

are the Green functions for propagation from chain i to 
chain j whereas G^± are the propagators for the sin- 
gle bands. Note that these relations hold only for the 



linearized model and therefore are valid only for small 
magnetic fields. Large magnetic fields of order 1/2 flux 
quantum per plaquette in the original lattice model affect 
the curvature of the dispersion relation and can eventu- 
ally lead to a gap!^ Such a field would be unrealistically 
large for a condensed matter realization of such a model; 
while this limit may be possible in certain superstruc- 
tures or optical lattices, this case will not be discussed 
here. 



B. Interactions 

We introduce a screened, short-ranged electron- 
electron interaction via two terms, one describing forward 
scattering 



H intjs = \^2 J dx 



(2.13) 

with the electron density n fJ/Cr (x) = ^^(x)^ ^ a (x), and 
one term describing backward scattering 

H int , bs = - J2 j dx^-^gr'^-^^' 

dx^l^-gf^-^^'- (2-14) 

fitra' 

Note that the interaction as shown here is similar to that 
of a single chain spinful Luttinger liquid. In consequence, 
the band index a can also be seen as a pseudospin; the 
unusual looking g\ term (corresponding to spin-flip scat- 
tering in the spinful Luttinger liquid analogy) appears 
as particle number in individual bands is not conserved 
by the interaction. Each of the coupling constants can 
be further split into components acting within the same 
band (a — a') and different bands (a = — a'). 

At low temperatures, the interplay between the dif- 
ferent interaction terms lead to a rich phase diagram 
of strong coupling phasesp^2l j n w hich the pseudospin 
excitations acquire a gap. However as temperature is 
raised above this gap the system simply behaves as a 
Luttinger liquid, exhibiting pseduospin-charge separation 
and power-law correlations. In this regime, the only role 
of the back-scattering terms g\ and g\ is to renormalize 
the forward scattering terms, and so they may be ne- 
glected. Furthermore, in this regime we will neglect the 
band dependence on the interaction terms and set 



gT 



gT 



9- 



(2.15) 



The reason for this is twofold - firstly, these interaction 
terms correspond to scattering events with low momen- 
tum transfer so in any realistic situation they are likely to 
be similar anyway; secondly and more importantly, the 
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most crucial role of the interaction terms and the con- 
sequent Luttinger liquid state for the current work will 
turn out to be the spin-charge separation, which already 
occurs in the single parameter model so extra complexity 
is not needed. 

It will be convenient to introduce the dimensionless 
interaction strength 



g/irvF- 



(2.16) 



As we will show, this one parameter is responsible for 
both the exponent of renormalization of disorder and de- 
phasing timep^ and will be assumed to be small a«l. 



C. Disorder 

The final part of the Hamiltonian is disorder scatter- 
ing Ha is = J2 a =i 2 / dxU a {x)n a {x) where the white-noise 
disorder potential is taken to be uncorrelated between the 
two chains and chosen from the usual Gaussian distribu- 
tion 



(U n {x)U n ,(x')) 



5(x — x')8 n 



2ttvt, 



(0) 



(2.17) 



Here, v = (ttvf) 1 is the density of states per chain, and 
T t r is the (bare) transport scattering time. The disorder 

is considered to be weak, e/r t ^ 3> 1. Since forward 
scattering does not affect transport, we only consider 
backward scattering off impurities, so U(k ~ 0) = Ut = 
and U(k ~ 2k p) — Uf,. Thus we can write the disorder- 
induced term in the Hamiltonian as 

H dis = J2 I dX ( U b^R.n^L,n + ' *I'L, <!' /,'. 

n=l,2 ^ 

( 2J8 ) 

The backscattering amplitudes are correlated as 

(U b , n (x)U b>n ,(x')) = and (t7 b ,„(*)£/*„, (a/)) = 

(U n (x)U n '(x')), furthermore we assume that the impu- 
rities scatter only intra-chain, but not inter-chain. 

At the most elementary level, the disorder gives rise 
to diffusive behavior of the system and one observes a 
Drude conductivity 



e 2 v F 

OD = 2 Ttr 



(2.19) 



however in the interacting syst em the scat tering time is 
renormalized by the interactio n 1 7 ! 13 ! 14 ! 27 ! 28 ! and therefore 
gains a temperature dependence 



r (0) I T 
tr 1 A 



(2.20) 



where a is just the dimensionless interaction strength 
introduced in Eq. (2.16), and T is temperature. 



In fact, for the present model, this equation is strictly 
only valid for the limit T > v^Afc i.e the temperature 



is greater than the splitting between the bands. This is 
because below this temperature, the renormalization of 
certain inter-band scattering processes will be cut-off by 
the band splitting and not temperature^ giving rise to a 
magnetoresistance effect even at the level of the Drude 
conductivity!^ Furthermore, this band dependence on 
renormalization will mean that the impurity scattering 
no longer remains diagonal in the chain basis. However, 
for weak interaction this renormalization is weak, and one 
can ignore this asymmetric renormalization so long as 
we stay away from the very low temperature T -C vpAk 
limit. In everything that follows, r tr is assumed to take its 
already renormalized, and therefore weakly temperature 
dependent, value. 

Having introduced the Hamiltonian, we now proceed 
to a calculation of the weak localization correction, which 
gives the dominant behavior of the magnetoresistance. 



III. CALCULATION OF WEAK 
LOCALIZATION CORRECTION 

In the absence of interaction, disorder leads to full 
localization of all electronic states in one- and two di- 
mensional systems; leading to zero conductivity at any 
temperature. However, this is no longer the case in the 
presence of interaction. Inelastic electron-electron scat- 
tering leads to a finite dephasing length 1$, which causes 
a finite weak localization correction to the conductivity, 
AawL- At high enough temperatures, the correction is 
small compared to the Drude conductivity ao, but at 
lower temperatures AawL becomes of the order of op>- 
At this point, the system crosses over to the strong lo- 
calization regime. 

The weak localization correction has previously 
been calcu lated for a single-channel spinless Luttinger 
liquidMi3! anc j a \ so f or the spinful caseP^ We follow a 
similar approach to calculate the weak localization cor- 
rection for the spinless disordered two-chain ladder which 
is layed down in Sec. [IT] We restrict the calculation to the 
limit of strong dephasing, t^/t <C 1, which is always the 
case for either sufficiently high temperatures or for suf- 
ficiently low disorder. In this limit, the disorder may be 
treated perturbatively, and we only have to retain the 
shortest possible Cooperon, that is, the one with 3 im- 
purity scattering lines 

The difference between the pure one-dimensional sys- 
tem and the two-chain ladder is that the electrons can 
perform loops which enclose a finite area, thus, weak lo- 
calization is sensitive to a magnetic field, which leads 
to magnetoresistance. An example of such a process is 
shown in Fig. [3]for a specific setup of impurity positions. 
Of course, at the end one should average over all possible 
impurity configurations. 

The weak localization correction depicted schemati- 
cally in Fig. [3] corresponds formally to calculating the 
diagrams in Fig. |4j the essential point being that the di- 
agrams should be fully dressed by interactions. We now 
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FIG. 3. Example process for scattering on 3 impurities where 
a finite area is enclosed. Stars indicate impurities, the red line 
shows a possible way an electron can take in order to scatter 
on all three impurities and the shaded region denotes the area 
which is enclosed by the loop. 




a) 



b) 



FIG. 4. Leading order diagrams to the weak localization cor- 
rection in the limit r^/r <C f . Dashed lines describe backscat- 
tering off impurities, the current vertices are dressed with dif- 
fusons and the solid lines are Green functions with disorder- 
induced self-energies. The diagrams are understood as fully 
dressed by electron-electron interactions. 



calculate this in two different ways - firstly in Sec. Ill A 
we treat the interactions phenomcnologically assuming 
that their only role is to introduce a dephasing time 
into the problem which is introduced by hand. The ad- 
vantage of this approach is transparency - one ends up 
with an analytic answer and can understand the different 
geometrical processes that give rise to different terms in 
this answer. Secondly in Sec. MB we calculate the di- 



agram within a full microscopic approach by employing 
the technique of functional bosonizationP^l By matching 
the phenomenological approach with the microscopic one 
in Sec. IIII CI we show that the former calculation is for 
the most part a very good approximation, and we find an 
explicit expression, Eq. ( |3.28 ), for the relation between 
the dephasing time t$ and the interaction strength ex. We 
also discuss an effect that goes beyond the phenomeno- 
logical picture and depends in a more essential way on 
the nature of the correlations in the Luttinger liquid - 
namely magnetoresistance oscillations. 



A. Phenomenological approach 

We start by calculating the weak localization correc- 
tion in a phenomenological way, closely following the 
method presented in Appendix B of Ref. [T3J The basic 
idea is that one considers the retarded Green function 
G.r/l.ct( w 7 k) to be of the non-interacting form (bu t with 
a renormalized scattering time, as discussed in Sec. II C I: 



Gr/l,c(u, k) = [w- v F (±k - k F ,a) + i/4rtx] 



(3.1) 



the numerical factor in l/4r tr comes from the specifics 
of the model with only backscattering, and expresses the 
total scattering rate in terms of the transport scattering 
rate. One then adds a phenomenological dephasing time 
introduced as an additional source of decay via the 
substitution 

Gr/l.<t{u, k) ->• [lj - v F (±k - k F ,a) + i/ 4r tr + i/^P 1 • 

. (3-2) 

The crucial difference between the transport lifetime r tr 
and the dephasing time is that the latter should be 
included only for Green functions running within the weak 
localization loopW^ In this sense, the inverse dephasing 
time which accounts for all inelastic processes acts as an 
effective lower cutoff on allowed momenta in the weak 
localization correction, as it shouldP^ 

The leading order diagrams for the weak localization 
correction in the limit t^/t 1 are shown in Fig. [4] The 
two diagrams [4Jd and [4j; sum up to the value of diagram 
[4^,, so the total correction is given by 



Aa#™ = 2<r C3 . (3.3) 
ac'3 is the fully crossed diagram [3£i and is given by 
de ( df € \ f dQ 



&C3 



-2 (ev F y 



2tt 



de 



2tt 



J{Q), (3.4) 



where the 2 comes from the two possibilities to choose the 
chiralities at the current vertices, the minus sign reflects 
the different signs of the two current vertices and f e is 
the Fermi distribution. The function J(Q) is given by 



J(Q) 



l 



E 

ijklm 



dp 

2n 



dpi 
2tt 



dp 2 
2tt 



x G 



^\p)G { t k \ Pl )Gf l \p 2 )Gf 



\-P + Q) 



Gf\p)Gr\-p 2 + Q) 



y(l,k) 



xG% J \-Pi + Q)G2' m, (-P + Q) 



OV 



(3.5) 



where we sum over all possible chain setups of the im- 



purities. From Eq. (|2.11|) we can transform the Green 
unctions 

J(Q) = 



functions into the band representations and find 
1 



(2irvT tI ) 3 



E 



pijkl 



dp f dpi 
2ir / 2tt 



ijkli'j'k'V 

x G R ^{p)G Lj {pi)G^ k {p 2 )G LJ (~p + Q) 

X G R ^(p)G L ,y(~p 2 + Q) 

x G R>k , (-pi + Q)G L ,i> (-p + Q) 



dp 2 
2ir 



(3.6) 



The sum over the different chains as well as the prefactors 



in Eq. (2.111 have been absorbed in P\^} k iy, the sum in 
Eq. (|3.6| is over the different band setups. The exact 

— ' ki n 

definition of P^j, k ,i, is given in Table I 
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(1 - b 2 ) 3 /8 


{i,j,k,i} 


{i, k,j,i} 


(l + 6 2 ) 3 /8 


{i, -i} 




(l + & 2 )(l-& 2 ) 2 /8 




{i,~3, 




-i, -i} 




(b 2 -fe 6 )/8 




{i,3, 




{i,3,i, -i} 


{i, -i} 


-b 2 (l-b 2 ) 2 /8 


{i,3,3,i} 


{i,~3, 




-i, -i} 


{i,3,3, -i} 


6 2 (l-6 2 ) 2 /8 


{i,j,k,-i} 


{i, k,j, -i} 




{i,3,3,i} 


{i,3, 






{i,-j,3,i} 



TABLE I. Values of the factors Pfg\ for a given set of indices 
{7} and {8}. All latin indices can exhibit the values + and 
— . Combinations not listed here result in Prgy = 0. The two 
sets {7} and {8} can be switched without changing the value 
of P {l} 



The current vertices should be dressed by the inser- 
tion of diffusons; for the backscattering only model this 
amounts simply to the replacement evp — > evp/2. Fur- 
thermore, the Green functions connected to the current 
vertices can be decomposed, 



2ir tT 



1 - 2iT tr v F [kF,a - k F ,<r'] 



(3.7) 



which leads to 

J(Q)= E WWW 

ijkli' j'k'l' 



(3.8) 



The prefactors of Eq. (|3.7|) have been absorbed in the 

rfijkl 
i' j'k'l' 



definition of P^-,u/i, and the functions F(Q) are given by 



1 



dp 
2^ 



G Pli (p)(G? w -(Q-p))* (3.9) 



aid F$ = {F%)\ 

Note that Eq. ( |3.8[ ) describes only the loop between the 
impurities, the diffusive motion has already been taken 
care of by Eq. (3.7 1. Consequently, it is at this stage we 



make the substitution ( |3.2[ ) and include the phenomono- 
logical dephasing time in the Green functions. Hence the 



F(Q) functions may be evaluated 
1 



Ftf(Q) = F^(Q) = 

F+ R {Q) = 
F^(Q) = 



2r tr l/2r tr + l/r^ + iv F Q' 



(3.10) 



1 



1 



2r tr l/27t r + l/r + iv F (Q - Ak) ' 

1 1 
2^ l/2r tr + l/r + iv F (Q + Ak) ' 

where Afc = hp _ — kp+. 



Combining (3.101, (3.8) and ( 3.4 ) we therefore find for 
the weak localization correction to the conductivity in 
the limit -C r tr to be 



A phcn _ _1 (l) ( 



K phcn (b, VpTfAk). (3.11) 



Here, — e 2 v F T tr /Tr is Drude conductivity of a sin- 
gle chain - this is introduced for convenience, the actual 
Drude conductivity of the two-leg ladder (2.19) is sim- 



ply given by <yp> — 2<7p)'- The dimensionless function 
-Kphen contains the interesting functional dependence of 
the weak localization correction and is given by 



-Kphen(MF7>Afc) 

1- - 



1 + {v F T^Ak) 2 (1 + ( VFTlP Ak) 2 ) 2 



-9 



15 



-9 

1 + {vpT^Aky 

32 

4 + (vpT^Ak) 2 
15 



1 



{v F T<j>Ak) 2 
-48 



-2 

(1 + (v F ^Ak) 2 ) 2 
192 

(4 + (vpr^Ak) 2 ) 2 
-2 

(1 + (vpT^Ak) 2 ) 2 
-256 



4+(w F r Afc) 2 (4 + (vpT^Ak) 2 ) 2 



-7 



+ 



1 + {vpr^Ak) 2 (1 + (v FT ^Ak) 2 ) 2 

16 64 

+ 



4 + {vpr^Ak) 2 (A + ivpr^Ak) 2 ) 2 



(3.12) 



with b defined in Eq. (2.8). The properties of this result 
will be examined in detail in Sec. |IV[ but first we will 
compare this with a direct microscopic calculation. This 
will allow us both to relate the dephasing time introduced 
here to the microscopic parameters of the model as well 
as to establish limits on the validity of the above formula. 
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B. Microscopic approach 

Having treated the electron-electron interaction qual- 
itatively in the last section via a phenomenological de- 
phasing time, we now turn to a true microscopic cal- 
culation using the method of functional bosonization. 
This method was introduced in Refs. [35] and [37J for the 
case of a clean Luttinger liquid and further developed in 
Refs. [34j l38ti40l The framework has recently also been 
extended to the diso rdere d case, where the transport 
properties of a spinlesa^EH and spinfufS! disordered Lut- 
tinger liquid have been studied. Functional bosonization 
has the advantage that it preserves both fermionic and 
bosonic degrees of freedom. This property is extremely 
helpful in the case of disordered Luttinger liquids, since 
interaction and disorder can be treated on an equal foot- 
ing. Here, we extend the previous work on single (but 
possibly spinful) channel Luttinger liquids to the case of 
a spinless two-leg ladder. This extension is simplified by 
remembering that the two chains may be thought of as 
pseudo-spins, which allows us closely to follow the for- 
malism developed in Ref. Q3J 

The basic principles of functional bosonization as ap- 
plied to disorder diagrams are summarized in Appendix 
|A} In essence, one factorizes the full Green function 
G M (x, t) in real space into a non- interacting part (x, r) 
and an exponent of an interaction (bosonic) correlator, 



Gp{x, t) = gl(x, t) exp r)] 



(3.13) 



Here, /x is the chirality, and the interaction correlator 



Bfj,fj.(x,T) which is given by Eq. (A6). The crucial point 



now is that in our interaction model Eq. (2.151, only for 



ward scattering terms with zero-momentum transfer are 
retained. Hence all information about the original band 
structure of the ladder model - that is, the splitting of 



the bands by the transverse hopping and the external 
magnetic field - is encoded in the free Green function 
<7°(x, r), leaving the interaction propagators unaffected 
by the geometric details. The splitting of the two bands 
in the two-leg ladder (giving each band its own Fermi 
momentum fci? j(T ) is accounted for by an additional phase 
factor in the free Green function: 



9°,ai X > 



(3.14) 



where g^(x,T) is the free Green function of the single 
channel Luttinger liquid, given by Eq. ( A5 ) in Appendix 
[A] The propagators in the chain basis are then given 
by the usual linear combinations, Eq. (2.11). Hence the 



only difference between the present situation of a two- 
leg ladder (in a magnetic flux) and the single chain with 
spin is the addition of linear combinations of these phase 
factors. This property allows for an easy adaption of the 
methods in Ref. EJ we therefore outline the calculation 
referring to the original paper for full details. 

When calculating diagrams within the functional 
bosonization approach, the relevant interaction propaga- 
tors must be added between every pair of vertices (see 
Appendix |A 2\ . For observable quantities (which are 
given by closed fermionic loops), this general principle 
amounts to retaining factors of 

Q(x,t) = exp [B rr (x,t) + B ll {x,t) - 2B rl (x,t)] . 

The Q(x, t) functions stem from the disorder, every 
pair of backscattering vertices at points (xn,tn) and 
(xjsf , tjv' ) is accounted by a factor Q(x, r) (with x — 
xn ~ %N' , t — t~n — tjv' ) if the chiralities of the incident 
electrons are the same, and by a factor of Q~ 1 (x,t) if 
they are different. 

Using these relations, the leading order weak localiza- 
tion correction in the functional bosonization scheme is 
given by 



^wl° = ^l<ev F ) 2 ]T P$ l m ,[ dndfidrid^dndn J d Xl dx 2 dx 3 

\ m ijkli' j'k'V 

x \g%(xi - x 3 ,ri - f 3 )Q~ 1 (xi - x 3 ,n - f 3 )][g° Lj (x 2 - x\,t 2 - t 1 )Q~ 1 (x 2 - x\,t 2 - t%)] 
x i9°R,k( x 3 ~ x 2,t 3 - t 2 )Q~ 1 (x 3 - x 2 ,t 3 - T 2 )][g° L (x 1 - x 3 ,fi - t 3 )Q~ 1 {xi - x 3 ,fi - T 3 )] 
x [dR,k'( x 2 - xi,f 2 - f 1 )Q~ 1 (x 2 - xi,f 2 - n)][gl tj ,(x 3 - x 2 ,f 3 - f 2 )Q- 1 (x 3 - x 2 ,f 3 - f 2 )] 
x Q(xi - x 3i tx - t 3 )Q(xi - x 3 ,n - f 3 )Q(xx - x 2 ,ti - f 2 )Q(x 2 - xi,t 2 - n) 
x Q{x 3 - x 2i t 3 - f 2 )Q(x 2 - x 3 ,t 2 - f 3 )(3~ 1 (0, n - fi)Q _1 (0,r 2 - f 2 ) 

x Q-^O^-fa^^-xs^fa,^ \ , (3.16) 



r 



where the space-time points correspond to the ones in Fig. [5j the factor 4 = 2x2 comes from the 2 different 




are defined as 



dx z dT l e- iamT *e- [lx <*- x > l+lx e- x > l]/l : 



x 9°nJx a ~ Xi,T a - Ti)g° Jxi - x fj 



(ftira ■ at 



-Tfl) 

(3.17) 



and 



''(s.T-a.T^n™) = W^'^X,^,^,-^). (3.18) 



FIG. 5. Three-impurity Cooperon, crosses denote impurities, 
solid lines represent free Green functions. The space-time 
points of the impurities are denoted by N = (xn,t~n) and 
TV = (xN,Tfj)- Each pair of backscattering vertices has to be 
accounted by a factor of Q or Q _1 (see text). 



choices to choose the chirality at the current vertices and 
the summation over diagrams (see Fig. [4]) , L is the system 
size and the P^^^y are the same geometric factors as for 
the phenomenological case, given in Table |IJ The factor 
p/2Z stems from the impurity correlator defined in eq. 



(2.17), using the (rcnormalized) mean free path I — vpTt T 



The factors W. 



in/ out,H.[ L 
i,3 



(x, r, t', tt m ) come from the in- 



tegration over external coordinates and times of the two 
Green functions attached to the current vertices [com- 
pare to Eq. (3.7 1 in the phenomenological picture], and 



While it is not difficult to evaluate the full expression 
for W" l /° ut ' R/L (x,T,T',n m ) (see Ref. 14), it can be im- 
mediately simplified by noticing two things. Firstly, all 
summands in Eq. (3.16) where i ^ i' or I ^ V vanish 
(see Table so only the W™ / °"*" R/i (a;, t, t', fi m ) with 
the same band indices survive. Secondly, in the strong 
dephasing limit which we are considering, only the lead- 
ing order expansion in x/l is needed; longer loops being 
exponentially supressed. Finally, by including the vertex 
corrections (which amounts to replacement of the total 
scattering rate 2vp/l by the transport rate vf/1) we find 
thatP 

Wl"'^(x,T a ,Tp,Cl m ) = 

_ ic tiik FA x sgn(f2 m ) ^ e ~in m r a _ e ~in m r^ 



v F /l 



(3.19) 

From Eqs. (|3.14|) and (|3.19|) we find that the sum in 



Eq. (3.16) simply consists of phase factors e lkF - i ' x and 



the factors P\^} k 'V Performing the sum results in the 
function 



T WL (b, Akx a , Akx b , Akx c ) = ^{fe°[5 + cos(2Afcr a ) + cos(2Afoz; fc ) + cos(2Afcr c )] 

+b 2 [ - 9 + 6(cos(Afcr a ) + cos(Akx b ) + cos(Afca; c )) 

- l(cos(2Afca; Q ) + cos(2Afcx h ) + cos(2Afcr c )) 

- 2(cos(Afc(x a — x b j) + cos(Afc(x a + x c )) + cos(Ak(x b + x c )))] 
+b i [15 - 8(cos(A£;x a ) + cos(Afcx h ) -I- cos(Afca; c )) 

- l(cos(2Afca; Q ) + cos(2Afcx h ) + cos(2Afcz; c )) 

+ 4(cos(Afc(ir a — x^)) + cos(Afc(x a + x c )) + cos(Ak(x b + x c )))] 
+b 6 [ - 3 + (cos(2A£;x a ) + cos(2Afcz; 6 ) + cos(2Afez; c )) 

+ 2(cos(Afc(x c — x a )) + cos(Ak(x a + Xb)) + cos(Ak(x c — Xb))) 

- 2(cos(Afc(x a - x b )) + cos(Afc(x a + x c )) + cos{Ak(x b + x c )))]\, (3.20) 



r 



where Afc = kp _ — kp + is the splitting of the Fermi points, b is given in Eq. (2.8) and the new variables 

x a — X\ X3 , 
x b = x 3 - x 2 , 

x c = x\ — X2, (3-21) 
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which clearly satisfy x a + x b 



The remainder of Eq. (3.16) (i.e. everything which 
carries no band index) is identical to the expression for 
the weak localization correction of the spinful Luttinger 
liquid (apart from a factor of 2 from spin degeneracy 
which now appears in the geometrical factor Twl), which 
was solved in Ref. 14 1 The essence of the calculation 
(see Ref. Q3] for details) relies on noticing that in the 
weakly interacting limit a <C 1, most of the Q factors 
cancel, the only interaction correlators remaining being 
those corresponding to the Green functions themselves. 
Putting all of these factors together, we find that the 
weak localization correction is given by 



Aa W L _ y j_ v 



n^o\ n m 32Z 4 (|fi m | +v/iy 



E 



pl/T poo 

/ dT a f a dT b f b d,T c T c I dx a dx b dx c 6(x a + x b - x c ) 
Jo Jo 

x exp[z(Sl„ + e n )(r a + n + t c ) + ie n (f a + f b + f c )] 
xT WL (b, Akx a , Akx b , Akx c )G R {x ai T a )G L (x a ,f a ) 



x G R (x b , T b )G L (xb, n)G R (x c , t c )G l (x c , t c ) 



(3.22) 



where G p without band indices refer to the full interact- 
ing Green functions but without the fast oscillating e lkpX 
factors; and are given in Eq. ( A17 ) in Appendix A 3 This 



is conveniently written in terms of the Green functions 



in the mixed (cc, e) representation (see Appendix A 3 1, 
which gives: 



AawL ,. f 2nT v 4 
= hm 



o£) n^o\ n m 32Z 4 (|O m | +v/l) 2 

dx a dx b dx c S(x a + x b — x c ) 

x G r + (x a ,ie n + iCt m )G r + (x b ,ie n + itt m ) 
x G r + (x c , ie n + iil m )G°L (x a , ie n ) 
x G a _( )G a _(x c ,ie n ) 



x TwlQ>, Akx ai Akx b , Akx c ) 



in m -vf2+i0 

(3.23) 



where G r + and G°_ are given in Eq. (A19). 

Finally, performing the analytical continuation iCl m — > 
0, + iO, we find in the DC limit fi — 5-0 



A ^micro 
A<T WL 



where L 



1 



.(i) 



AkL 



, (3.24) 



vf/chT is the electron-electron scattering 



length and K micm (b, 7) is given by 



v in icro (Pi T) — 
7T 



dx 



dyx 



4 J^oo cosh^(7rz) J Jo 
x H(x, z)H(y, z)lZ(x + y + xy, z)x 
x T WL (6, 1 L{x), 1 L{y), 7 (L(x) + L(y))) . 

(3.25) 

Here L(x) = ln(l + x), 7 = Akl ee /2 and 

72,(35, z) =2^i(l/2 + iz, 1/2, 1; -x)x 

x 2 F 1 {l/2-iz,l/2,l;-x), (3.26) 

where 2 F\{a, 6, c; z) is the Gauss hypergeometric func- 
tion. 



The integral in Eq. (3.25) must be evaluated numeri- 



cally. In the next subsection, we will show this result, and 
compare it to the far simpler expression (3.12 ) calculated 
via the phenomenological approach. 



C. Comparison between phenomenological and 
microscopic calculations 

So far, we have calculated the weak localization correc- 
tion of the spinless two-leg ladder in a phenomenological 
and a microscopic approach, the results are given in Eqs. 
(3.11) and (3.24). Both results have a similar structure, 



even the functions A' p h cn /micro (h, 7) have the same struc- 
ture, 

K p {b, 7) = q p fi{i) + b 2 q p . 2 {-i) + b 4 q pA (~f) + b & q pfi {l), 

(3.27) 

where p = phen/micro, 7 = v^r^Afc in the phenomeno- 
logical and 7 = Akl ee /2 in the microscopic case and the 
Qp.i (l ) can be obtained from the different terms in Eqs. 
(3.12) and (3.25). Note that while this may look like the 



first few terms in a power series expansion in (dimension- 
less) magnetic field b, the series naturally terminates due 
to our limitation to loops containing only three impuri- 
ties. 

The two-parameter functions K p (b,^/) may thus be ex- 
pressed in terms of a linear combination of four func- 
tions of a single parameter, making comparison between 
the phenomenological and microscopic results easy. The 
only thing to fix is the relationship between the dephas- 
ing time introduced by hand into the phenomenolog- 
ical model, and the electron-electron scattering length 
l ee that was a natural parameter in the microscopic 
approach. In fact, this relationship may already be 
fixed by looking at the limiting value of AawL when 
Afc = 0, i.e. the limit of zero inter-chain hopping. Fixing 



0) = Aer^™(A£; = 0) yields the relation- 



ship 



T 4> 



1.09 



.L 



1.09 



v p iraT 



(3.28) 
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FIG. 6. [Color online] Comparison between gphcn.i^FT^Afc) (black solid line) and numerical values of g m icro,i(wFT"</>Afc) (red 
dots) as defined by Eq. (3.271. The blue curve is a spline fit to the numerical data points. The microscopic q's have been scaled 
so that qo = 2, in accordance with the phenomenological result. This scaling also gives the relationship between dephasing 
time T,), and electron scattering length l ee which fixes the horizontal scaling of the microscopic result (see main text). Panel (a) 
after appropriate scaling also has physical meaning as the weak localization correction to the conductivity of the ladder in the 
absence of an external magnetic field - see section |TV| of the main text. 



We can now plot the four functions functions 5 p hen,i('y) 
and <7micro,i(7) in Fig. [6j Note that having established 
the relationship between and l ee , this comparison has 
no fitting parameters — the numerical values obtained 



via Monte Carlo evaluations of the integrals (3.25) are 



simply plotted alongside the analytic curves obtained 
from p. 12 ). This comparison demonstrates two things 
- firstly, the overall shape of the curves is very well cap- 
tured by the phenomenological approach; and secondly, 
the microscopic calculation shows that these functions 
also contain an oscillating component not captured in 
the simple phenomenological calculation. 

The fact that these oscillations are not seen in the 
phenomenological picture means that they are related 
to features in the correlations of the two-leg ladder be- 
yond dephasing, and nothing to do with normal magneto- 
resistance oscillations in two dimensions which origi- 
nate from the Landau levels (in fact, the basically one- 
dimensional two-leg ladder does not exhibit Landau lev- 
els) . In Appendix |Bj we demonstrate that the impor- 
tant feature giving rise to the oscillations is that of 
pseudospin-charge separation, i.e. the fact that the two 
leg ladder has two (collective) modes propagating with 
different velocities. In the next section, we show how 
these oscillations translate into oscillations in the mag- 



nctoresistance, and give a more physical description of 
their origin. 



IV. DISCUSSION OF WEAK LOCALIZATION 
AND MAGNETO-CONDUCTANCE 

We turn now to an examination of the properties of 
the weak localization correction given by Eq. (3.24), 



tx,f), (4.1) 



where as usual, a^y is the Drude conductivity of a single 
chain and K{T^,t^_, f) is either the scaled microscopic 
result K = (1.09) _2 X m i cro or the phenomenological ap- 
proximation to this K = ifp hcn- These functions are 
given by Eqs . ( 3.12) and (3.25 1, where b and Afc are given 
by Eqs. ( |2.8| and (2.10); and for the microscopic case the 
relationship between and l ee is given in Eq. (3.28 1 . Ul- 



timately, K(r^,,tj_, f) depends only on the dimensionless 
ratios 



#"(t0,*x, f) = K[ T^tx, 



t± 



(4.2) 
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however to keep the discussion more transparent, we re- 
tain the more physical original parameters. 

The weak localization correction of the spinless Lut- 
tinger liquid is given bjG^ 



LL 



(4.3) 



which differs from Eq. (4.1) only by the function 
K{r^,tx, /)• In order to investigate the properties of 
the weak localization correction, we examine the func- 
tion K(t^, t±, /) and relate the results to the case of a 
single chain. We discuss first the smooth part of the 
weak localization correction, for which we can use the 
analytic results from the phenomenological calculation 
-Kphen) which matches the microscopic value in all limit- 
ing cases. We then discuss the oscillations on top of this, 
for which we must switch to the results of the numerical 
integral, K miclo . 



A. Weak localization correction in the absence of a 
magnetic field 

In the absence of interchain hopping wc find 

K phcD {T^,t ± =0,/) = 2, (4.4) 

which means that the weak localization correction of the 
two-leg ladder is twice as large as the correction of the 
single chain. This is the result one would expect in the 
case of two decoupled chains. Furthermore, this result 
is independent of magnetic field which must be so as for 
two decoupled chains, the magnetic field is a pure gauge. 
However, this result is a non-trivial check on the calcu- 
lation, as for the case t± = in the band picture, zero 
magnetic field corresponds to b = 0, Afc = while finite 
magnetic field corresponds to b = 1, Ak — > oo. From 
Eq. (3.12), we then see that this gauge independence cor- 
responds to the statement 5 + 1 + 2 = 5 — 9 + 15 — 3. 
A simple check confirms that the related equality is also 
satisfied by K m i Cio . 

When interchain hopping is turned on but the mag- 
netic field is absent, we find 



K 



phcn 



(t 0) *x,/ = O) 



1 + (2^r ) 2 (1 + {2t ± T^y 



(4.5) 



which is smaller than 2 for non-zero t\T^. Hence, 
the weak localization correction of the two-leg ladder is 
smaller than the correction of two uncoupled chains. This 
result can be attributed to dimensionality, since localiza- 
tion effects become less pronounced with increasing di- 
mensionality of the system.^ In the limit of strong inter- 
chain hopping, T^t±_ > 1, we find A' phcn (r0, tj_, / = 0) -> 
5/4. The weak localization correction is still stronger 
than in the case of a single chain, but we cannot give a 
simple physical explanation of the factor 5/4. 



(2) 
(1) 



VfT± 



FIG. 7. Hopping from chain 2 to chain 1 in the continuous 
ladder model. Hopping can occur on a length scale of l±, as 
shown by the shaded lines. 



The weak localization correction in the absence of a 



magnetic field (4.5 1 coincides with the function we previ- 
ously called go j an d is plotted in Fig. [6^, as a function of 
VFT(j,Ak = 2T0ij^ when / = 0; which shows the smooth 
interpolation between the limits 2 and 5/4. The figure 
also shows the equivalent result from the microscopic the- 
ory, which more or less follows the smooth phenomeno- 
logical curve with some oscillations on top of this which 
die away as r^t^ — > oo. One can understand this with the 
following simple picture: the parameter r$ti_ is roughly 
speaking the number of times an electron can hop from 
one chain to the other before losing it's coherence. When 
this number is small, the fact that the charge and pseudo- 
spin modes propagate at specific different velocities gives 
rise to coherent oscillations in the weak localization cor- 
rection. If this number is too large, the average over 
different disorder positions means such oscillations are 
washed out and only the smooth part (as captured by 
the phenomenological theory) remains. 



B. Magneto-conductivity — smooth part 

We first consider the case of a strong magnetic field, 
where we find 



Kphen(T<t>,t±,f ->■ OO) ->• 2, 



(4.6) 



which is the same as the limit of two uncoupled chains. 
Surprisingly, the weak localization correction is stronger 
in the presence of a magnetic field than without it, so the 
magnetic field enhances localization. 

This rather counterintuitive result is a consequence of 
the specific geometric setup of the two-leg ladder, and 
can actually be understood rather easily, in a similar 
vein to the dying away of oscillations discussed above. 
A single electron propagating along one of the chains 
can hop to the other one on a typical length scale of 
l± = vpTj_, where we have introduced the average hop- 
ping time t_l = t^ 1 (see Fig. [7]). Since tj_ acquires a 
position-dependent Peierls phase in the presence of a 
magnetic field, the hopping parameter feels a different 
phase, depending where the hopping takes place. These 
hopping events interfere destructively and lead to an ef- 
fective suppression of interchain hopping, and hence is 
equivalent to the limit of two decoupled chains. This be- 
havior is already seen in Eq. (2.11 ), where the inter-chain 
Green functions vanish for large magnetic fields, b — > 1. 
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FIG. 8. [Color online] Plots of dimensionless magneto-conductivity &cj{t^, t±, f) against the generalized magnetic field, vpTrf /t± 
for (a) t±Tj, = 0.2, (b) t±T^ = 0.4 and (c) tj_T^ = 1.5. The black solid lines are obtained from the phenomenological calculation, 
while the blue dashed lines correspond to the equivalent microscopic results, approximating the numerical functions by the 
spline fits as obtained in Fig. [5] The black curves show the evolution from a monotonic to non-monotonic magneto-conductivity 
as t±T$ is increased, as discussed in Sec. |IVB| The blue curves show the magneto-conductivity oscillations on top of this, as 
discussed in Sec. IIV CI 



We now define the dimensionless magneto-conductivity 



Ao-wz,(T0,f j_,/) - Aa WL (T (j> ,t±,0) 
Aa WL (T <j) ,t ± ,0) 



(4-7) 

to study the behaviour of the weak localization correction 
for finite magnetic fields at a given T$t±. We have plotted 
5a(T t f 1 ,t±, /) for different values of ij_r^ in Fig. [8j 

For small values of t±r^,, the limit of two decoupled 
chains is achieved in a monotonic way. For large values 
of tj_T0, however, the behaviour is non-monotonic. This 
comes about from the competition of the above peculiar 
effect of magnetic-field enhanced localization, and the 
more usual decohering effect of the magnetic field that 
reduces quantum corrections to conductivity in higher 
dimensions.^. When t±T t f ) is large, the dominant con- 
tribution to the weak localization correction comes from 
loops in which the electron hops coherently many times 
between the two legs of the ladder - and hence for small 
fields such loops will feel the normal decoherence effect 
of a flux penetrating such loops. Only at larger magnetic 
fields does the peculiar enhancement of weak localization 
in this geometry take over through the effective decou- 
pling of the chains. 

A numerical analysis shows that the transition from 
monotonic to non-monotonic behaviour occurs in the 
phenomenological model at t±T<f, ~ 0.266. However, 
Fig.[8]clearly also shows that at these small values of tj_T<p 
the oscillations in magneto-conductivity arising from cor- 
relation effects are large; so the neat change in behavior 
predicted by the phenomenological calculation is not pos- 
sible to see. The oscillations will be discussed below. 



logical calculation - thereby demonstrating that these 
oscillations go beyond the single particle (with dephas- 
ing) picture, and depend in a more essential way on 
the nature of the strong correlations ubiquitous in one- 
dimensional interacting systems. In Appendix [B] we 
showed on a technical level that these oscillations have 
their origin in the pseudospin-charge separation; while in 
Fig. [8] we saw how this translates into oscillations in the 
magneto-conductivity which are particularly pronounced 
when i_i_T0 < 1. 

In fact, there have been a number of predictions of the 
effects of spin-charge separation on transport in clean 
Luttinger liquids. Some examples include time-resolved 
measurements of the propagation of pulsesp^l or interfer- 
ence effects in finite ringg^HSl w ith an Aharonov-Bohm 
flux. In fact, for sufficiently short ij_T0, the weak local- 
ization in our ladder system can be thought of as qualita- 
tively similar to a finite ring as in this limit one expects 
simple loops with only two inter-chain hopping events 
before the electrons become dephased. Therefore, the 
oscillations we predict in magneto-conductivity have es- 
sentially the same physical origin as the manifestation 
of spin-charge separation in Aharonov-Bohm ringsP^HlH 
It is an intriguing result that physically observable ef- 
fects of the interaction-induced correlations in the two 
leg-ladder system appear more naturally in transport in 
the disordered case than in the clean case. 

However, as t\jr$ becomes larger, many more coher- 
ent inter-chain hops of the electron take place, and the 
beautiful interference effects arising from two different 
velocities become washed out as previously explained. 



C. Magneto-conductivity — oscillating part 

We previously showed in Fig. [6] that the microscopic 
calculation yields oscillations in the weak localization cor- 
rection which were not seen in the simple phenomeno- 



D. Experimentally measurable conductivity 

In an experimental setup, the parameters easily varied 
are temperature and magnetic field, while one measures 
the full conductivity including the Drude component 



a{T, B) = a D (T) + Aa WL {T, B). 



(4.8) 
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FIG. 9. [Color online] Plot of full conductivity (including 
Drude component) as a function of magnetic field for various 
temperatures, from top to bottom T — 50, 40, 35, 32, 30K. 
The other parameters used in the plot are a = 0.01, A = 
10000Jf, v F = 10 8 cm/s, a y = lnm, t± = IK and l/r t ( r 0) = 
IK . The normalization of the vertical axis (Jd,o is the Drude 
conductivity in the high temperature limit. 



We remember that the weak temperature dependence of 
the Drude conductivity <jd comes from the Luttinger 
liquid renormalization of the transport scattering rate, 
while the strongest temperature dependence of the WL 
correction is through the temperature dependence of the 
dephasing time oc 1 /T, which not only enters the func- 
tional dependence on magnetic field but appears as a 
prefactor in the strength of the quantum correction. In a 
typical single experiment, other material constants such 
as interaction strength, or inter-chain hopping, will re- 
main constant. 

In order to relate to potential real experiments, we 
plot in Fig. [9] the conductivity against magnetic field for 
various temperatures, being careful to choose a param- 
eter range which is both physically reasonable, and in 
which the approximations in our calculations are valid 
- the exact parameters chosen are shown in the caption 
to the figure. As expected, the weak localization correc- 
tions (and hence the magnetic field dependence) become 
stronger as temperature is reduced. Furthermore, we see 
clearly the fact that, as a general trend, conductivity de- 
creases with increasing magnetic field - the effect we call 
magnetic field enhanced localization. We also see oscil- 
lations in conductivity as a function of magnetic field - 
which we remember are not due to conventional effects 
such as mesoscopic fluctuations (we average over disor- 
der in our calculations) or Landau levels (which do not 
exist in the two- leg ladder); but in fact have their origin 
in correlations, namely pseduo-spin charge separation. 

It is instructive to compare our results to the experi- 
ment on nanotubes by Man and MorpurgopS One should 
not compare details too closely as the geometric set-up 
of a nano-tube differs from that of a two-leg ladder; how- 
ever it is already striking how similar some of the fea- 
tures of our Fig. [9] and the experimental curve (Fig. 3 
of Ref. ITS]) look. It remains work for the future to adapt 
this calculation to the nanotube geometry relevant for 



this particular experiment; and to see whether this type 
of orbital magneto-resistance is applicable to single wall 
carbon nanotubes. 



V. SUMMARY 

In summary, we have calculated the weak localiza- 
tion correction to charge transport and the magneto- 
resistance for a model system of a disordered two-leg lad- 
der. We have limited ourselves to the high-temperature 
regime defined by the two conditions: 

1. Temperature T should be greater than any interac- 
tion induced gaps. 

2. Dephasing time oc (aT)^ 1 should be much less 
than the mean free time r tr of diffusion in the 
(weakly) disordered system. 

In this regime, we find that the weak localization cor- 
rection is a universal function of t\T^ and f/tx- Our 
results can be summarized as follows: 

1. In the absence of a magnetic flux, the weak local- 
ization correction to the Drude conductivity is a 
decreasing function of the size of the inter-chain 
hopping tj^T0 from 2 at t± — to 5/4 at t± —> 
00 in units of the single chain correction. This 
crossover is simply an aspect of the phenomena 
that as dimension is increased (in this case going 
from one pure one-dimensional chains to two cou- 
pled chains) , weak localization corrections decrease 
- this is plotted in Fig. [6^,. 

2. There are two main competing effects contributing 
to magnetoresistance. The first is the decoherence 
of weak localization loops by the magnetic field, 
leading to diminishing weak localization correction 
as flux is increased. The other effect is the ten- 
dency of the magnetic field to decouple the two 
chains, and thus enhancing the weak localization 
correction. For tj_T0 not too small, the two effects 
compete and lead to a non-monotonic magnetore- 
sistance - this is shown in Fig. [Hj 

3. When t±r^, is small, there is a third observable 
effect seen in the magnetoresistance, which is co- 
herent oscillations due to the different velocities of 
charge and pseudo-spin excitations in the two-leg 
ladder - this is also shown in Fig. [5J 

Finally, Fig. [9] puts all of these results together, and 
shows a typical example of what we predict would be 
seen in an experimental measurement on such a system. 

We expect that these results, though derived for a spe- 
cific toy model, are more general than this as we can 
identify the basic physical origin of each effect. In par- 
ticular, results 1 and 2 are even seen in a simple phe- 
nomenological model where the important features are a 
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phenomenological dephasing time and the system geom- 
etry. Our microscopic calculation backed up these gen- 
eral observations, as well as showing the beautiful ad- 
ditional effect of magnetoresistance oscillations arising 
from (pseudo) spin-charge separation, which we believe 
goes beyond the specific geometry of the two-leg ladder 
and may even be seen in e.g. nanotubes. 

We acknowledge useful discussions with A. Yashenkin 
and M. Schiitt, and support from the German-Israeli 
foundation. IVG acknowledges illuminating discussions 
with D. Polyakov during the preparation of Ref.H"4l about 
a related ladder model that turns out not to exhibit or- 
bital magnetoresistance. STC acknowledges early dis- 
courses with N. Hussey and A. Narduzzo which in part 
motivated the present work. MPS acknowledges useful 
discussions with T. Sproll. 



Appendix A: Functional bosonization 

In this appendix we present the method of func- 
tional bosonization for a disordered spinful Lut tinger liq- 
uid. The method has been used in the spinless^ESl and 
spinful 3 case; a more detailed derivation is presented in 
the given references. 



1. Basics 

The key ingredient of functional bosonization is a 
Hubbard-Stratonovich transformation which decouples 
the four-fermion interaction by introducing a bosonic 
field tp(x, r). This field has the correlation function 



{<p(x,T)<p(0,0)) = V(x,r), 



(Al) 



where V(x, r) is the dynamically screened interaction. 
This leads to a quadratic action for the Fermions; in- 
tegrating them out then gives the usual bosonized ac- 
tion - however this suffers from the disadvantage that 
fermionic operators are non-linear (and in some repre- 
sentations even non-local) in terms of the bosons, which 
creates difficulties in creating a diagrammatic expansion 
for disorder. In contrast, functional bosonization retains 
both the fermionic and bosonic degrees of freedom in the 
action; we will now outline how this overcomes the above 
difficulty and one can treat further perturbations to the 
system such as disorder. 

In one dimension, fermionic and bosonic fields can be 
decoupled in the action by the gauge transformation 

Vv( x > T ) = Vv( x ' T ) ex P[ i6 W x ' T )] j ( A2 ) 

where fj, = R/L is the chirality, a is the (pseudo-)spin 
and the phase 0^(x,t) obeys 

(d T -ifj,VFd x )6 fl (x, t) = lp^{x,t). (A3) 

When we write down the full Green functions in terms 



the Gaussian averages over the 6(x, r)-fields, we find [see 
e.g. Ref.|34] 

Gfj,(x, t) = g^x, t) exp [~B^(x, t)] . (A4) 

Here, the free Green function (ignoring for now phase 
factors e ±lkpX - see main text) are given by 



g^(x,r) = -fj,- 



iT 



2vf sinh [ttT(x/vf + *M r )] 
while the bosonic correlator B flv (x,T) defined as 

B^{x,t) = ([^(0,0)-^(x,t)]^(0,0)). 
is given by 

V(q,iQ n ) 



(A5) 
(A6) 



(A7) 



{v F q - iQ n )(±v F q - iQ n ) ' 

The remaining chiral combinations are given by 

B ll (x,t) = B rr (-x,t), B lr (x,t) = B rl (x,t). 

(A8) 

In this expression, f2„ = 2imT is the bosonic Matsub- 
ara frequency, and the dynamically screened interaction 
V(q,zQ n ) is given by 

„2 „2 U Q2 

% (A9) 



V(q,iil n )=g V lf. 



u 2 q 2 



where the charge velocity u is given by 



= Vps/l + 2g/nv F . 



(A10) 



Within our original model with a single interaction pa- 
rameter, g, the (pseudo-)spin velocity remains unrenor- 
malized and so is equal to tij?. 



Substituting Eq. ( |A9[ ) into Eq. ( A7 ) and executing the 
integral, one finds 



1 C^b 

Brr(x,t) = ~-\nr](x,T) - — ln?(x,r), 



Brl(x,t) 

where 



ln?(a;,r), 



(tvT/ A) 2 



(All) 



sinh [nT(x/u + it)] sinh [ttT(x/u — ir)] 
vp sinh [ttT{x/vf + it)] 



^ ' ' u sinh [ttT(x/u + it)] 
The prefactors a^ r in Eq. ( |ATT| > are given by 
(u- v F ) 2 



Ctb = 



a 
~2 



(A12) 



a, (A13) 



2uv F 2 2uv F 

where a — g/irvF is the dimensionless interaction 
strength introduced in Sec. |IIB| In the limit of weak in- 
teractions, the prefactors follow the hierarchy 



a b <C a r <C 1, 



(A14) 



of the fermionic operators given in Eq. ( A2 ) and perform 



since aj, is quadratic in the interaction strength, whereas 
a r is linear. 
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-ie fl (l) 



,;e B (2) 



-Bm( x 2-xi,n-T\) 




b) 



FIG. 10. Green function of a right-mover propagating from 
space-time point 1 = (xi,Ti) to 2 = {x2,Ti) before (a) and 
after (b) averaging over the 9(x, r)-fields. The solid line is 
the free Green function, the wavy line stands for the factors 
exp[— i9u(l)] and exp[i#R(2)]. After averaging over the 9- 
fields, the points 1 and 2 become connected, the wavy line in 



(b) stands for the factor e 



-S1.T2— n) 



2. Diagrammatic expansion 



Equation ( A2 1 suggests the following diagrammatic 



technique: The v(x, r)-fields get attached to the starting 
and ending points of the free propagator, drawn as wavy 



lines in Fig. 10 1. After averaging over the ^-fields, the 
two lines become co nnected and represent the correlator 
e~ B w - see Fig. 10 d. This diagrammatic technique can 
be extended to scattering off impurities in a straightfor- 
ward way. A backscattering process is shown in Fig.fTTfc, 
interaction is accounted by the 0-fields attached to the 
free Green function. Averaging results in connecting all 
wavy lines pairwise in every different way. 

In closed fermionic loops, the -B^-correlators appear 



only in the combination (see Fig. 11 ) 



M(x,t)=B rr (x,t) + B ll (x,t)-2B rl (x,t). (A15) 

As a consequence, each pair of backscattering vertices at 
space-time points N and N' contributes a factor of 

Q(x N - x n >,t n - tjv') = exp[M(x N - x n >,t n - t n >)] 

(A16) 

if the chiralities of the incident electrons at the vertices 
are the same, and a factor of Q^ 1 if they are different. 



3. Green function in (a;, e) representation 



Combing Eqs. (A4), (A5), (All I and (A12), and as- 



suming o-b *C 1 (weakly interacting limit) so that the <r 
term may be neglected, the interacting Green function 
takes the form 



G m (x,t) ps -fj, 



ttT 



ttT 



smh[nT(x/vF + «M T )] sinh[7rT(a;/u + i/ir)] 



(A17) 

1/2 



Within this approximation, the only singularity of 
as considered as a function of r is a branch cut between 
t = \x\/u and |x|/«i? in the upper of lower- half plane 



R 



e i[e fl (2)-e £ (2)] 

\ , 

— fex^ 

2 

(a) 



„iS L (3) 




(b) 



FIG. 11. (a) Backscattering process of a right-mover on 
an impurity (denotes by a cross), both Green functions are 
dressed by ^-fields, (b) The possible connections of the 9- 
fields for averaging in the simplest loop; which yields a Q 
factor, Eq. (A16l, associated with the pair of vertices. In a 



more complicated loop, each pair of vertices gives a factor Q 
or Q _1 (see main text). 



depending on the sign of fix. Choosing carefully the 
contour^ one can therefore represent the Fourier trans- 
form of the Green function with respect to r as 

l/T 

dre^G^x, t) = G^x, ie n ) - Gftx, ie n ), (A18) 

where e n = 2ix{n + V)T is the fermionic Matsubara fre- 
quency and 

G^x^en) = 0(e n )9(nx)g(\x\,\e n \) 

G«(x,ie n ) = 6{-e n )0{-nx)g{\x\, \e n \). (A19) 

Executing the integral gives G(x > 0,e„ > 0) asP 

1 e n 1 



^—e n x/ u-\-x/2l ee 

G(x,e n ) = . 2-Fi 

iJuvf 



2ttT' 2 



1: 1 - e 2x ' 1 ' 



(A20) 

where 2-P1 (a, b, c; z) is the hypergeometric function and 
lee — vp/aT is the electron-electron scattering length. 

For a full discussion of the properties of this Green 
function, see Ref. Q3] 
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Appendix B: Oscillations arising from different spin 
and charge velocities 

In this appendix, we show mathematically the relation- 



ship between the phenomenological calculation Eq. (3.11 ) 
and the microscopic result Eq. (3.24). By making this 



comparison, we isolate the source of the oscillations in the 
weak localization correction as seen in the microscopic 
calculation, and show that these may be attributed to 
pscudo-spin-charge separation. 

Starting from the integral for the microscopic calcu- 
lation, Eq. (3.251, we see that the 1/ cosh 2 (7T2:) factor 
means that the dominant contribution comes around 
z sa 0. Therefore using the expansion of the hyperge- 
ometric function 

2*1 (1/2, 1/2, 1; -x) « ln(16x) + o( X -^ 2 ) (Bl) 

and ignoring the logarithmic corrections to the power 
law, we see that this integral may be approximated as^ 



Aay/L 



dx a 



b -^Jl 
-e >ee x 



x T WL (Akx a , Akx b , Ak(x a + x b )) 

2 



K phen (b, Akl ee ), 



(B2) 



where we have returned to the original variables x a and 
xi, which are physically distances between impurities. We 
may now consider what will happen when the logarith- 
mic factors present in Eq. ( Bl I are taken into account. 



While these factors will change the true asymptotic be- 
havior of the Green function, this is actually irrelevant 
for evaluating the integral above as the contribution from 
these regions is exponentially small. We will shortly show 
that the slowly varying logarithm at intermediate length 
scales can actually be very well approximated by a slight 
modification of the power law in Eq. ( |B1[ ); which trans- 
lates in the integral above to a slight modification of l ee . 
In fact, fitting this power law at this intermediate length 
scale (see the in set t o Fig. 13 1 gives us exactly the mod- 
ification in Eq. (B2| of 



lee — ► 'd> — VpTA 



1.09L 



(B3) 



in full agreement with Eq. ( 3.28 ) in the main text. In this 



way, we have rederived the phenomenological result out 
of the microscopic model. It is worth noting that this 
single number governs both the vertical and horizontal 
scaling in all four graphs in Fig. [6| 

This may be seen even clearer by noticing that the 
above expansion essentially corresponds to writing the 
Green function of the electron as 



G(x a >0,e n «0)^e-^/ 2 ^ 



(B4) 



which is exactly the form used for the phenomenological 
Green function, Eq. (3.2 1. 




FIG. 12. [Color online] The differences between the micro- 
scopic and phenomenological results Aqi as a function of 
7 = VFT,pAk. Universality of the oscillations is seen when 
each Aq is scaled appropriately (see text). Specifically, the 
plot shows Ago (7/2) (blue solid), AffcM/o (red dashed), 
— Aq4(7)/8 (black dotted) and Aqe("/)/2 (green dash dotted). 



We now turn our attention to the oscillations, seen 



in the numerical evaluation of the integral (3.25) which 
clearly go beyond the approximations used above. We 
isolate the oscillations by looking at the difference be- 
tween the full calculation and the approximation above 
(or in other words, the phenomenological result), 



Aft 



<Zphc 



(B5) 



where i — 0,2,4,6. The oscillations (as a function of 
Aklee) in the final result come from integrals of the co- 



sine terms within the factor T\yl given in Eq. (3.20) 



Noting that there is a certain symmetry in T\yl between 
interchange of x a , x^ and x c , and singling out the simplest 
oscillating term for each power of b, we find that we would 
expect the dominant contributions to come from the term 
with cos(2Afcx a ) in q ; 6cos(Afcx a ) in q 2 ; — 8 cos(Akx a ) 
in g 4 ; and 2 cos(Afc(x c — Xf,)) = 2cos(A/ca; a ) in q e . Hence 
by scaling each Aqi by an appropriate factor, we expect 
to see that the oscillations in each q function are identical 
- which is demonstrated in Fig. [12] The slight differences 
follow from the sub- leading terms in T\yl- 

From this demonstration of the universality of the 
weak localization oscillations, we see that to capture their 
essence, we can suffice by approximating 

TwL(Akx a , Akxb, Ak(x a + Xb)) ~ cos(Afca; a ). (B6) 

It is then clear why the phenomenological approximation 
does not produce any weak localization oscillations, as 



dx a 



cos(Akx a ) 



1 



21 1 + (Ak lee/2) 



(B7) 



and the only role of the band-splitting is to renormalize 
the effective dephasing time. In fact, carefully collect- 
ing all such factors explains immediately the form of the 
phenomenological result, Eq. (3.12). Another important 



point can be drawn from the above formula however - 
which is that such contributions to the weak localization 
correction mathematically take the form of the Fourier 
transform of the single particle Green function (actually, 
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jtTx 



FIG. 13. [Color online] Comparison between the phenomeno- 
logical product of Green function, J- p hcn(x a ) (red dashed 
curve) and the true microscopic curve J-micro(xa) (black solid 



curve), as defined by Eq. (B8| The blue dotted curve shows 
the difference between them (multiplied by a factor of 5 for 
clarity). We see that as compared to the phenomenological 
dephasing model, the true microscopic model has enhanced 
spectral weight for distances of the order of a few times the 
dephasing length. The inset shows the comparison between 
the same microscopic and phenomenological functions on a 
longer length scale and a logarithmic axis. 



a product of six such Green functions, but this distinction 
is unimportant in the following qualitative analysis). To 
have oscillations as a function of Afc in the final result, 
the Green function should therefore have a peak at some 
finite value of x a /l ee . In order to see how this comes 
about, we consider the microscopic kernel Eq. (3.25) as 
a function of x at fixed y 



0, and compare this with the 
approximation above that yields the exponential - more 
precisely, we consider the following two functions: 



J~ micro (^a) ^ 



_>,. /,. / d n(e 2 ^-l,zfTZ(0,z) 



cosh (ttz) 



phen(^a) ^ 



_ p -2x a /l<, 



(B8) 



where TZ(x, z) is given by Eq. (3.26 1. The value of used 
in J-phen is chosen to best match the exponentially decay- 
ing tail of .Fmicro at intermediate length scales (see the 



inset to Fig. 13 ), and matches the 1$ as fitted in the main 
text. We emphasize that this is however a phenomeno- 
logical fit to intermediate length scales which dominate 
the integrals for the weak localization correction; and not 
the true asymptotic of the microscopic function, which 
involves power law prefactors. 

We now turn to shorter length scales, of the order of 
several l ee ; in this range, the two functions are plotted in 
Fig. [T3j along with a magnified view of the difference be- 
tween them. We see that as compared to the phenomeno- 
logical model, the true microscopic model has correla- 
tions which give an enhanced spectral weight at distances 
of order a few times the dephasing length. Having iden- 
tified mathematically this feature that is responsible for 
the weak localization oscillations, we can finally ask the 
question: which property of the correlations in the inter- 
acting two leg ladder are responsible for this enhanced 
spectral weight? 



FIG. 14. [Color online] Real time propagation of spectral 
function |Ss[<5 p hcn/micro(£, t)]\ f° r interaction strength a = 0.2 
and times nTt — 0.1,0.75,1.5,3,5,7.5 (delta functions are 
broadened to show the correct spectral weight under them). 
The microscopic model (solid lines) propagates with two 
peaks due to the two velocities, but the spectral function re- 
mains large between these peaks, explaining the extra spec- 
tral weight at intermediate distances unaccounted for by the 
phenomenological model (dashed lines). 



To answer this question, we consider the real time 
propagation of the true microscopic Green function, 
which is given by an analytic continuation of Eq. ( A17) 
and yields 



^micro (*£; ^) 



(ttT)* 



1/2 



(1 + a) sinh7rT(x - t) sinh%T( T ^ t - t) 



(B9) 



and compare it to the real time propagation of the phe- 
nomenological Green function, given by a (finite temper- 
ature) Fourier transform of Eq. (3.2), 



<?phen(M) OC 



nT e 



v/2 



sinhTrT^-t) 



(BIO) 



In both cases, we have set Up = 1, and have made use 
that both the inelastic scattering length and the differ- 
ence in pseudo-spin and charge velocities are proportional 
to the interaction strength a. 

The time evolution of these functions is shown in 
Fig. [14] In the microscopic model, it is easy to see the 
propagation via two wavefronts, reflecting the different 
charge and pseudo-spin velocities; however it also has the 
curious feature that there is significant spectral weight 
between the two peaks which is not accounted for in the 
phenomenological picture. It is exactly this feature which 
after integrating over time gives enhanced spectral weight 
at intermediate distances, and is therefore responsible for 
the oscillations in the weak localization. At longer dis- 
tances x a l ee , this physical property is also responsible 
for the power law enhancement of J^nicro as compared to 
■^phen _ but as these distances contribute an exponen- 
tially small amount to the weak localization, we do not 
comment any more upon this asymptotic behavior. 
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